clc; close('all'); clear;
load('Datos11164');
T = 0.1;
tn = 50;        % Segundos
N = tn / T;     % Numero de muestras para tn Segundos
time = a_ent_esc(:,1);
in = a_ent_esc(:,2);
out = c_pert_esc(:,2);
txtC = {'Perturbacion','Respuesta en C','Recta',...
    'Estimacion del Modelo G2(s)*G3(s)'};
[Tu, Ta, K] = identificar(time, in, out, N, txtC);
T1 = Ta / 2.99;
T2 = 1.2*T1;
s = tf('s');
G2G3 = K /((1+T1*s)*(1+T2*s));
